Split transition in ferromagnetic superconductors 
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The split superconducting transition of up-spin and down-spin electrons on the background of 
ferromagnetism is studied within the framework of a recent model that describes the coexistence 
of ferromagnetism and superconductivity induced by magnetic fluctuations. It is shown that one 
generically expects the two transitions to be close to one another. This conclusion is discussed in 
relation to experimental results on URhGe. It is also shown that the magnetic Goldstone modes 
acquire an interesting structure in the superconducting phase, which can be used as an experimental 
tool to probe the origin of the superconductivity. 
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I. INTRODUCTION 

Recently, the coexistence of ferromagnetism and super- 
conductivity has been observed in a number of materials, 
including, UGe 2 r*£ URhGe^ and ZrZn 2 i^ The experi- 
ments so far have ascertained the presence of bulk super- 
conductivity, and various thermodynamic and transport 
properties have been measured, but little is known yet 
about the detailed nature of the superconducting state. 
An obvious possibility is spin-triplet pairing induced by 
ferromagnetic fluctuations, but other mechanisms have 
also been proposed^ One interesting aspects of the ex- 
periments, which may provide a clue about the origin and 
the nature of the pairing, is that the superconductivity 
is observed only on the ferromagnetic side of the mag- 
netic phase boundary. This is in sharp contrast to early 
theories of superconductivity induced by ferromagnetic 
fluctuations, which predicted that this type of supercon- 
ductivity would be equally strong on the paramagnetic 
side of the magnetic phase boundary^ However, in a re- 
cent paper^ to be referred to as I, it has been shown 
that the presence of magnons in the ferromagnetic phase 
can lead to a drastic enhancement of the superconducting 
T c compared to the paramagnetic phase. Clearly, more 
properties of the superconducting state must be studied 
in order to differentiate between different possible pair- 
ing mechanisms. In this context, an interesting result is 
the recent measurement of the specific heat in URhGe, 3 
Taken at face value, this experiment shows a single transi- 
tion into a superconducting state, and a low temperature 
specific heat that is linear in the temperature, which sug- 
gests that a fraction of the electrons remain in a Fermi- 
liquid state at low temperatures. Similar behavior has 
been observed in UGe 2 £ 

To date all of the theoretical work has assumed only 
the simplest type of superconducting order for a given 
pairing mechanism, with the main goal being to deter- 
mine the phase boundary for superconductivity. For ex- 



ample, in I the present authors assumed an ordering of 
spin-triplet Cooper pairs with spins oriented in the direc- 
tion of the magnetization. Let us denote the gap function 
for this ordering by At. Previous work on Helium-3 in a 
magnetic field 9 suggests that at some point, the Cooper 
pairs with spins oriented opposite to the magnetization, 
characterized by a gap function A^, will form as well. 
More generally, the complete phase diagram for these 
systems will likely be complicated and involve ferromag- 
netism coexisting with several types of superconducting 
order ii& 

Our goal in this paper is three-fold. First, we develop a 
formal theory that enables us to consistently describe fer- 
romagnetism coexisting with two types of superconduct- 
ing order. The superconductivity in our theory is caused 
by ferromagnetic fluctuations. Our goal is to derive an 
equation of state, analogous to the strong-coupling or 
Eliashberg equations of conventional superconductivity, 
that describes both up-spin and down-spin superconduct- 
ing order, as well as a consistent magnetic equation of 
state in the presence of superconductivity. A major com- 
plication is that, because the superconductivity is itself 
caused by a fluctuation effect, a simple mean-field the- 
ory is not sufficient, and fluctuations need to be taken 
into account. The fluctuations that cause superconduc- 
tivity within our theory are described by the spin suscep- 
tibility tensor. Our second goal is therefore to develop 
a theory for the spin susceptibility, both in a pure fer- 
romagnetic phase, and in coexisting ferromagnetic and 
superconducting phases. The pairing potential for su- 
perconductivity involving only At and A^ is given by 
the longitudinal susceptibility, xl- However, the trans- 
verse susceptibilities, xt, are also needed because they 
enter the normal self energies, and because they couple 
to xl via mode-mode-coupling effects. In fact, as shown 
in I, within the framework of our theory it is this modc- 
mode-coupling mechanism, which exists only in a mag- 
netically ordered state, that causes the superconductivity 
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FIG. 1: Schematic phase diagram in the temperature- 
pressure plane showing a paramagnetic phase (PM), a nor- 
mal conducting ferromagnetic one (NCFM), a superconduct- 
ing ferromagnetic phase with up-spin pairing only (SCFM I), 
and a superconducting ferromagnetic phase with both up-spin 
and down-spin pairing (SCFM II). The two superconducting 
phases have been drawn to end in the same point for simplic- 
ity, but this will not necessarily be the case. 



to be observable only in the ferromagnetic phase. We will 
see that the magnons, which are described by XTj have 
an interesting dispersion relation in the superconducting 
state and, under certain conditions, can become effec- 
tively massive. A striking feature of the mass is that it is 
proportional to the inverse square of the magnetization, 
and thus singular for small magnetization. We give esti- 
mates showing that this result should be observable with 
neutron scattering. 

Our third goal is to discuss the phase diagram for su- 
perconductivity and magnetism, based on the above re- 
sults. In a theory that allows for pairing of both up-spin 
and down-spin electrons, one expects a phase diagram 
that is qualitatively shown in Fig. ^ With decreasing 
temperature, one oberves a transition from a paramag- 
netic (PM) phase to a normal conducting ferromagnetic 
(NCFM) one. With decreasing temperature, the sys- 
tem enters a ferromagnetic superconducting state where 
only up-spin electrons are paired (SCFM I). Finally, with 
further decreasing temperature, down-spin electrons are 
paired as well (SCFM II). In general the coupled mag- 
netic and superconducting equations of state are very 
complicated to solve. For the points we want to make, 
however, a complete solution is not needed. Our aim here 
is to compute, for our proposed pairing mechanism, the 
relative magnitudes of the transition temperatures T c f 
and T c | for the phase boundaries between the NCFM 
phase and the SCFM I phase, and between the SCFM I 
phase and the SCFM II phase, respectively. Generically 
we find that these transition temperatures are close to 
one another. Experimentally this suggests that, for ex- 
ample, any specific heat measurement should find two 
closely spaced transition signatures. This result is in 



conflict with the naive interpretation of the specific heat 
measurement in URhGe noted above. A possible alterna- 
tive interpretation of the data consistent with our theory 
will be discussed below. 

The plan of this paper is as follows. In Section ITT1 we 
develop a formalism that allows for a consistent descrip- 
tion of all components of spin-triplet superconductivity, 
induced by ferromagnetic fluctuations, in the presence of 
long-range ferromagnetic order. In Section IlIII we calcu- 
late the magnetic susceptibility in the ferromagnetic state 
in the presence of a non-unitary superconducting order 
parameter. In Section IIVI we solve the strong-coupling 
equations for superconductivity and determine the phase 
diagram containing phases with pure ferromagnetic or- 
der, ferromagnetic plus spin- up superconducting order, 
and ferromagnetic plus both spin-up and spin-down su- 
perconducting order, respectively. We discuss our results 
and their experimental implications in Section^] In Ap- 
pendix we augment our microscopic approach with a 
more general Landau-Ginzburg- Wilson theory and dis- 
cuss the soft-mode structure of the magnetic supercon- 
ducting phase. In Appendix [5] we relate the physical 
spin susceptibility to the magnetization fluctuations that 
occur most naturally in the theory^ Some of our results 
have been reported before in Ref. [jji 

II. A FIELD-THEORETIC APPROACH TO 
SUPERCONDUCTIVITY AND MAGNETISM 

A. The model 

Our starting point is the same model of interacting 
electrons as in I. That is, we consider free electrons with 
a static, point-like spin-triplet interaction with ampli- 
tude r t . For simplicity we ignore the spin-singlet in- 
teraction. We do not include an explicit Cooper chan- 
nel interaction; the pairing interaction will be generated 
by magnetic fluctuations. We then decouple the spin- 
triplet interaction by means of a Hubbard-Stratonovich 
transformation. Bilinear products of fermionic fields are 
constrained to a bosonic field Q by means of a Lagrange 
multiplier field A with a transposed A T , and we inte- 
grate out the fermions. The resulting action is given by 
Eq. (2.12a) of I, with T s = 0, 

A[M,g,A] = i Tr h^G^ 1 + \/%Tt~f ■ M — A T ) 

+Tr(A£) ~ J dx M(x) ■ M{x) . (2.1) 

Here Gq 1 is the bare inverse Green operator, 

G^ = -a T + 7 o(V 2 /2m e + M) , (2.2) 

where /i is the chemical potential and m c is the elec- 
tron mass. We denote the real-space and imaginary- 
time coordinates by x and r, respectively, combine these 
into a four- vector x = (x,t), and use the notation 
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j dx = J v dx J^ T dr, with V and T the system's vol- 
ume and temperature, respectively. M is the Hubbard- 
Stratonovich field whose expectation value is propor- 
tional to the magnetization m, and 7 denotes three com- 
ponents of a four-vector of 4 x 4 matrices, 

(70 = 7) = (°~3 ® cr , o 3 <g> o-i,(To O 0-2,0-3 <g> a 3 ) , (2.3) 

with 01,2,3 an d o the Pauli matrices and the 2x2 unit 
matrix, respectively. G(x, y) is a 4 x 4 matrix field which, 
in contrast to M, depends on two space-time variables. 
Its expectation value determines the various Green func- 
tions of the electron system. Specifically, 

{Gn{x,y)) = -{Qas{x,y)) = \ {^{x)ip^{y))^ 



1 



G r (x-y) 



(2.4a) 



(Q22(x,y)) 



-{Ga{x,y)) = - W>1 (a#| (?/))</> 



G^x-y) 



(2.4b) 



are the normal Green functions, while 



{Qia(x,y)) 
(02i(x,y)) 
(G3i{x,y)) 
(G42(x,y)) 



1 



(^(x)^(y))^ = -F+(x-y) , (2.4c) 



WdxWiMh = o f j ( x ~y) '( 2 - 4d ) 



(^(x)^(y))^ = -F r (x-y) , (2.4e) 
(TPi(x)i>i(y)h = lFi(x-y) , (2.4f) 



are the anomalous ones. Here 



and 



denote 



averages with respect to the effective action A and the 
underlying fermionic action, respectively. Since we do 
not allow for mixed pairing of up- and down-spins, which 
one expects to be strongly suppressed, these are the only 
nonzero Green functions. We decompose M and G into 
their expectation values and fluctuations 



M(x) = my/Tt/2 z + SM(x) 
G{x,y) = (G(x,y)) +5G(x,y) 



(2.5a) 
(2.5b) 



Here m is the magnetization, which we assume to be in 
^-direction. For the Lagrange multiplier field A we write, 
in analogy to Eq. (|2.5b|l . 

A(x,y) = \(x - y) + 6A(x,y) . (2.5c) 

The matrix elements of A represent the normal self ener- 
gies, 

Ml{x-y) = -A 33 (a; - y) = E T (a; - y) , (2.6a) 
^22{x-y) = -Xu(x-y) = T,i(x-y) , (2.6b) 

and the anomalous ones, 



Xl3{x 


-v) 


= A^(x 


-y) , 


(2.6c) 


X 2 a(x 


-y) 


= A+(x 


-y) , 


(2.6d) 


hi(x 


-y) 


= A T (x- 


-y) , 


(2.6e) 


\ i2 {x 


-y) 


= A l(^- 


-y) ■ 


(2.6f) 



Despite the formal similarities in our treatments of M 
and G on one hand, and A on the other, it is important 
to note that A is not the expectation value of A. Rather, 
it will be determined self-consistently by the method ex- 
plained below. 



B. Expansion in powers of the fluctuations 

We now expand the action in powers of the fluctuations 
8M and 5G, as well as the quantity <5A. To this end it is 
useful to define an inverse Green operator 



G- 1 =Go 1 +mT a3 -\ T 
The zeroth-order or mean-field action then reads 



(2.7) 



1. 



. >T m^ + Tr(A(g» + -TrlnG- 1 . (2.8) 



This mean-field action or Landau theory, and its gener- 
alization to a Landau-Ginzburg- Wilson (LGW) theory, 
are interesting in their own right, and we discuss them 
in Appendix 1X1 To bilinear order in the fluctuations, it 
is obvious from Eq. (|2.1|l that 6A couples to both 5G 
and SM . This coupling can be eliminated by first shift- 
ing SM, and then 5A. The diagonalized Gaussian action 
then takes the form 



A {2) = - {SM\x~ l \SM) - ^ (<5A|r|<JA) 
+l(5G T \T-i\SG T ) . 



(2.9) 



Here we employ an obvious scalar product notation that 
implies summation and integration over all discrete and 
continuous indices, respectively, of the various fields, x" 1 
is an inverse Gaussian magnetic susceptibility defined by 



where 



Xil{x-y) = 5 ij 5(x-y) + T tX %{x-y) , (2.10a) 

X%{x -V) = \* (G(y - x) 7i G{x - y)j 3 ) . (2.10b) 

In Appendix \B\ we show that x is directly proportional 
to the magnetic spin susceptibility X s m & Gaussian ap- 
proximation. The matrix T, which determines both the 
6 A and the SG propagators, is given by 



r = r (D +r (2) 5 



(2.11a) 



where 



r a p, a 'f3'( x y^ x 'y') = ^G apl {x~y')G alfs {x' -y) , 

(2.11b) 

x (7i) 7< 5 G S p(z-y)xtj(z- z')G a > Y {x' - z') 



x (7i)y ff ' G s , ,{z' -y') 



(2.11c) 
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We now consider the non-Gaussian terms in the ac- 
tion, which also are affected by the shifts that diagonalizc 
A^ . The terms linear in the fluctuations are represented 
graphically in Fig. [3 Explicitly, we find 



FIG. 2: Contributions to the linear part of the action, ■ 
Solid lines denote 8M, dashed lines 5A, and dashed-dotted 
lines 8Q. 



Notice that the 5Q vertex in Eq. (|2.9(l is minus the inverse 
of the 8 A vertex. This property will be important in what 
follows. 



J 



A {1) = (a | 8M) + (b | 6 A) + (c | 8Q T ) , (2.12a) 



with 



and 



'2T t 



T 

-mS l3 + — Tr ( Gh/i 



(2.12b) 



b a p{x,y) = {G T ) a p{x - y) - ^G al3 (x - y) + V2T\ J dx'dy' ^ ^ G aj (x - x') ( 7i ) 75 G S p(x' - y) Xij (x' - y') a 

ij -y 6 

c a p{x,y) = \ af 3(x-y)+ / dx'dy 1 ' ^2 b a > p (x' ' ,y')T~} pl af3 (x'y' \xy) . 



(2.12c) 
(2.12d) 



For the cubic terms wc find 



A {3) = J dxdydzJ2a ljk (x,y,z)8M i {x)SM J (y)SM k (z) + J dx dx' dx" dy" ^ ^ j3'^ ij ' (x, x'; x" , y") SM^x) 

ijk " ij a/3 

xSMj(x') 5Gjp{x", y") + (other terms) , (2.13a) 

The vertices read 



and 



where 



a ijk (x, V,z) = - (2r t )( 3 / 2 ) tr [g(z - x) 7l G{x - y) 7j G(y - z) lk 



(3'^; 3 (x, x'; x", y") = dz dz' ]T a ,^ (x, x'; zz') T^ ap (zz' , x"y") 

J a>0' 



(2.13b) 



(2.13c) 



f3' a p 3 {x,x';x" ,y") = /3 Q/3 U (x, x'\ x", y") + | y/T\ J dz dz' ^ ^ o^y (x, x', z) {'fi') > a , Xi'j'( z ' ~ z)G all3 {z' - y") 

i'j' ot'0' 

xG a p,{x" - z') , (2.13d) 



with 



f3 a p ij (x, x'- x", y") = -r t [G(x" - x) 7l G(x - x 1 ) 7j G(x' - y") 

I 



a/3 



(2.13e) 



We note that (3 is the original (SM) 2 SA vertex. Via the shifts of SM and SA that diagonalize A^ it generates 
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FIG. 5: 
order. 



The superconducting equation of state to one-loop 



FIG. 3: Contributions to the cubic part of the action, _4' 3 '. 
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FIG. 4: The magnetic equation of state to one-loop order. 

the (5M) 2 8g vertex (3". The 'other terms' in Eq. I|2.13aj) 
we will not need explicitly. Their structure is displayed 
in Fig. 

We are now in a position to determine the magnetic 
and superconducting equations of state. In I this was 
done by means of integrating out the fluctuations and 
minimizing the free energy with respect to the order pa- 
rameters. This procedure is not straightforward, 12 and it 
becomes more confusing the more order-parameter fields 
one needs to consider. For our present purposes, where 
we have one magnetic and two superconducting order 
parameters, we therefore prefer a generalization of Ma's 
methodic 



C. The magnetic equation of state 

Following Ma^ we determine the magnetic equation 
of state by requiring 



(5M(x)) = 



(2.14) 



This is a formal expression for the exact equation of state. 
By expanding the action in powers of the fluctuations, as 
we have done above, it can be evaluated order by order 
in a loop expansion. To one-loop order, we obtain the di- 
agrams shown in Fig.0] Consider the last two diagrams. 
The original action, Eq. JOJ, contained no 8Q legs at 
cubic order; these were produced by the shift 



sa -> r^sg 1 



(2.15) 



that decouples SA and 8Q. As a result, the vertices in 
these two diagrams, which in Fig. 0] are denoted by v 



and /x, respectively, are multiplicatively related by two 
factors of T -1 . Symbolically, 

v = nT- 2 . 

Together with the relation between the 5 A and 8Q propa- 
gators that was mentioned after Eq. (|2. 1 lc|) , this implies 
that the last two diagrams cancel each other: 

5MvY~ l -8MfiT = . 

Clearly, this mechanism is not restricted to these partic- 
ular diagrams, but yields the following general diagram 
rulei^ 

Rule: 8A loops and 8Q loops cancel each other. 

This is the reason why we did not need to explicitly de- 
termine the 'other terms' in Eq. I|2.13a|) . An evaluation 
of the first two diagrams in Fig. 0] yields 

mS i3 = ^Tr (jiG^ + y J dx dy^Xjkix - y) 



jk 



tr 



7iG(-x)7jG(x - yh k G(y) 



(2.16) 



The first term on the right-hand side represents the 
mean-field magnetic equation of state, which was derived 
(by a different method) and discussed in Iji^ while the 
second term represents one-loop fluctuation corrections. 



D. The superconducting equation of state 

We now determine the superconducting equation of 
state by requiring 



(8G T (x,y)) =0 



(2.17) 



Taking into account the diagram rule from the preceding 
subsection, this condition is graphically represented in 
Fig. El A calculation yields 



Xotp{x-y) = - J dx'dy'^2 

a' P' 



8 a 'p'(x' - y') 



- 1 - J dx"dy"Y,Xi j {x'-y')0 i i, p ,{x'y',x"y" 

ij 

^K'0',a i x "y"^y) ■ ( 2 - 18a ) 
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with x from Eq. (|2.10a|l . (3 from Eq. (|2.13e|l . T from Eqs. 
(UTTTJ, and 



(2.18b) 



Here we have used the magnetic equation of state, Eq. 
(I2~m 

The equations (|2.18|) for A must be supplemented by 
a relation between G and (£7). We stress again that 
A 7^ (A), so one must not require (SA) = O. 16 - Rather, 
we calculate (G) directly. Going back to the underlying 
fermionic formulation of the action, it is easy to show 
that 



(G T )(x,y) = ±(G M ,A(^y)) 



(2.19) 



where Gm,a is G, Eq. I|2.7|l . but with m and A replaced 
by the full fluctuating fields M and A. In evaluating 
(Gm.a{x, y)}, we restrict ourselves to one-loop order, and 
also to linear order in the magnetic susceptibility x, as 
contributions quadratic in x are indistinguishable from 
two- loop contributions. We find 



{9 T )= l -G + 0{x 2 ,2-\oo V ) 



(2.20) 



The first term in brackets in Eq. H2.18a|) thus vanishes, 
and the second term we again evaluate to linear order 
in x- Wc finally obtain the superconducting equation of 
state in the form 



Kxp(x - y) = r t y^Xijiv- 

— L J 0a 

(2.21) 

For later reference, we note the following. We are an- 
alyzing the equation of state in a loop expansion, and 
our treatment constitutes a systematic one-loop evalu- 
ation. However, since the superconductivity does not 
occur at all unless one goes to one-loop order, while the 
magnetism appears already at zero-loop order, the result- 
ing equations of state are not on the same level physically. 
Specifically, the magnetic equation of state, Eq. I|2.16|l . 
contains fluctuation effects, while the superconducting 
one, Eq. I|2.21|l . does not, except for the magnetic fluctu- 
ations that cause the superconductivity in the first place. 
In fact, the generalized Eliashberg equations that follow 
from Eq. H2.21fl are analogous to conventional Eliashberg 
theory, which neglects all superconducting fluctuations. 



E. The Eliashberg equations 

Writing Eq. I|2.21|l explicitly yields the desired Eliash- 
berg equations for superconductivity induced by mag- 
netic fluctuations. By using Eqs. I|2.7|l and l|2.3|) we ob- 
tain a set of coupled equations for the matrix elements 



of A, Eqs. i2~Hl) . 

A T (fc) = T t f XL(«-*)A T (g)/d T (g) , (2.22a) 

Jq 

A x (fc) = r t [ XL (<Z-fc)Ai(9)M(<7) , (2.22b) 

Jq 

£ T (A0 = r t / XL(.q-k)Gj\q)/di(q) 

Jq 

+2T t f [XT,+ (g-fc)+iXT,-(?-fc)] Gl 1 (q)/d l (q) , 
Jq 

(2.22c) 

Si(fc) = r t / XL(g-fc)G- 1 ( <Z )/d i (g) 

Jq 

+2T t f [XT,+ (9-fc)+iXT,-(«-fc)] Gj\q)/d 1 (q) . 
Jq 

(2.22d) 

The A + obey the same equation as the A. Here we have 
performed a Fourier transform to fermionic Matsubara 
frequencies oj n = 2irT(n + 1/2) and wave vectors k, and 
we use the notation k = (iuj n ,k), J k = T^2 n (l/V)^2 k . 
We have introduced 

d a (k)=G~ 1 (k)G- 1 (-k) + A a (k)At(k) , (2.23) 

and the GZ 1 , are the inverse 'normal' Green functions 



(see Eqs. (|23U|) . $F7\. and (|2~Ial l2~4bl . 

G~ 1 {k) = iuj n - Cfc.cr - S<r(fe) 



(2.24) 



with £fc a = k 2 /2m c — fx — aS. Here a =t, J.= ±, and 5 = 
Tt m is the Stoner gap or exchange splitting. Finally, we 
have used the fact that the magnetic susceptibility tensor 
X in the presence of a magnetization has the structural^ 



XT,+ XT- 
-XT- XT.+ 

o o XL/ 



(2.25) 



This structure holds in general, and in particular for 
the explicit approximate expression for x given by 
Eqs. (|2.1U|) . In a superconducting phase, with non- 
vanishing gap functions, x will depend on the gap. This 
gives rise to a complicated feedback mechanism that 
is characteristic of any purely electronic mechanism for 
superconductivity^ We will discuss this feedback in the 
following sections. 



III. THE MAGNETIC SUSCEPTIBILITY IN 
THE PRESENCE OF SUPERCONDUCTIVITY 

The Eliashberg equations (|2.22|) require the mag- 
netic susceptibility x as input, just like the Eliashberg 
equations for conventional superconductivity require the 
phonon propagator as input. There are two possible at- 
titudes one can take at this point. In principle, one could 
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use experimental results for Xi m analogy to experimen- 
tal phonon spectra being used as input for solving the 
conventional Eliashberg equations. Sufficiently detailed 
information for \, however, is not available. Alterna- 
tively, one can calculate or model x m an effort to con- 
struct a self-contained theory. This is what we will do 
now. For the purpose of determining the phase diagram 
for up-spin and down-spin superconductivity, we need x 
in two phases. The up-spin superconducting T c is de- 
termined by x m the normal conducting ferromagnetic 
phase, NCFM in Fig. ^ This was discussed in I, and we 
briefly recall the result in Sec. IIII Al For the down-spin 
superconducting T c , we need x m the phase that has both 
magnetic and up-spin superconducting order, SCFM I in 
Fig. H This is discussed in Sec. IIII Bl 



effect that couples xl to xt and leads to a longitudinal 
susceptibility that diverges at k — everywhere in the 
ferromagnetic phase. 20 As was discussed in I, the one- 
loop expression 

Xh\ k ) = ~^Xh{k) / [XT,+ 0- q)XT,+ (q) 

+ XT,-(*-g)XT,-(9)] XL (fc) , (3.2b) 

takes this effect adequately into account. The one-loop 
approximation 

XL(k) = X ( L ) (k)+x[ 1 \k) , (3.2c) 

thus correctly reflects the behavior of xl at small wave 
numbers and frequencies. 



A. Normal conducting ferromagnetic phase 



The Gaussian theory, Sec. IIIBI yields an explicit ex- 
pression for Xi y iz., Eqs. Q2. !()() . For the normal conduct- 
ing ferromagnetic phase, this was evaluated in I. For the 
transverse susceptibility tensor at small wave vectors and 
frequencies, the result is 



XT,+(k) 



1 



5/Ae F 

l-t \iQ, n /Ae F + (5/2e F )b T {k/2k F ) 2 
1 



Xt,- (k) 



in n /4e F - (S/2e F )bT(k/2k F y 



1 



(3.1a) 



-iS/4e F 

1 - t V^n/4e F + (<V2e F )&T(k/2fc F ) 2 
1 



+ - 



n n /4e F ~ (<5/2e F )& T (k/2fc F ) 2 



(3.1b) 



Here k = (ifl n , k) with f2„ = 2-kTti a bosonic Matsubara 
frequency, t = 1 — 2N F Tt is the mean-field distance from 
the magnetic critical point, with N F the density of states 
per spin at the Fermi surface, and e F is the Fermi energy. 
6t = 1/3 in our free-electron approximation, but more 
generally it is a number of order unity. This result dis- 
plays the magnons, or magnetic Goldstone modes, that 
are a consequence of the spontaneously broken spin ro- 
tation symmetry in a ferromagnetic phase. It provides 
a qualitatively correct expression for the transverse spin 
susceptibility in such a phase. 

For the longitudinal susceptibility, the Gaussian ap- 
proximation yields 



XL(k,zO) = 



l-t 



a h \t\+b L (k/2k F y 



(3.2a) 



with cil and 6l constants of 0(1). Popular model cal- 
culations give ol = 5/4 and b L = 1/3. 19 In contrast to 
the transverse channel, however, Eq. (|3.2a|l is not quali- 
tatively correct. The reason is the mode- mode coupling 



B. Superconducting ferromagnetic phase 

We now consider the magnetic susceptibility x i n the 
ferromagnetic phase with Aj =^ 0, = 0, SCFM I in 
Fig. ^ We first consider the Gaussian approximation 
defined by Eqs. H2.10|) . and discuss the validity of this 
approximation later. In terms of Green functions, the 
five nonzero matrix elements of xZ- read, 



X^{k) = 



1+ 2 Ft 



Gn(q + k) Gn(g) 



+G 22 (q + k) G 22 (q) - G 13 (q + k) G 3 i(q) 

+(*;-> -fc)} , (3.3a) 
1 + T t / G u (q) \G 22 {q + k) + G 22 (q - k) 

Jq 

Xn\k) , (3.3b) 



G U (q) G 22 (q + k)-G 22 (q-k) 



(3.3c) 



We need the susceptibility only at zero external fre- 
quency, and we perform the integrals in an approxi- 
mation that neglects the normal self energy as well as 
the frequency dependence of Af. That is, we approx- 
imate A-f(fc) ~ Af k z , with k z the z-component of the 
unit wave vector. It is convenient to do the summa- 
tion over frequencies first, and to replace the wave num- 
ber summation by an integral over £ q . The calculations 
are very similar to those of susceptibilities in an s-wave 
superconductor^ 

One readily finds that, as in the case of s-wave super- 
conductors, X33 — 0) is independent of the gap. 
For xl — I/X33 1 m Gaussian approximation, Eq. (|3.2a|) 
is therefore still valid. In the transverse channel the sit- 
uation is more complicated. It is illustrative to first con- 
sider the case of zero external frequency and wave num- 
ber. From Eq. (|3.3b(l . and using the magnetic equation 
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of state, Eq. in zero-loop approximation (i.e., ne- We see that, within the framework of our approxima- 



glecting the second term on the right-hand side) , we have 



(3.4) 



X22 (k = 0) = 1 + 2I\ / G 11 (q)G 22 (q) 



|A T (<z)| 2 



tions, the transverse susceptibility has a mass propor- 
tional to A| . For small values of the frequency, the wave 
number, and Af , we thus obtain, instead of Eqs. (|3.1|) . 



S/4e F 

1-t \in/4e F + {S/2e F )b T [(k/2k F ) 2 + (A T /2c5) 2 g(5/2T)} 
-iS/4e F ( 1 
1-t \in/4e F + (5/2e F ) b T [(k/2k F ) 2 + (A T /25) 2 g(S/2T)} 

I 



(8^-5) 
+ (S^ -5) 



(3.5a) 
(3.5b) 



Here 



g{8/2T) 



-6<5 



^f|A t P Jq 
In the limit T <^ e F we obtain 



a^Hg^)! 2 ^) 



(3.5c) 



IV. THE STRUCTURE OF THE 
SUPERCONDUCTING TRANSITION 

We now are in a position where we can solve the gap 
equations (|2.22|) for the down-spin superconducting crit- 
ical temperature, T c f . We will do so in an approximation 
that is analogous to the one employed in I for T c f . 



j(y) = dx 
Jo 



1 



x(l — X 2 ) 
where a = 0.8525 . . .. 



tanh(xy) = 



In y if y 1 

ay 



2 if y <C 1 



(3.6) 



The above results have been obtained in an approx- 
imation that neglects all fluctuations of the supercon- 
ducting order parameter, see the remark at the end of 
Sec. Ill D)l . This raises the question whether the mass is 
generic, or a result of our approximations. To answer this 
question we consider, in Appendix IA 21 a LGW theory 
that treats magnetic and superconducting fluctuations 
on equal footing. Both a general symmetry analysis and 
an explicit calculation within the framework of the LGW 
theory show that the result obtained above is correct for 
all wave numbers only in the limit ^a/^t™ — > oo, where 
£a and £ m is the superconducting and magnetic (zero- 
temperature) coherence length, respectively. For a finite 
value of this ratio, xt as a function of the wave number 
displays a shoulder, but still diverges for asymptotically 
small values of |fc|. We conclude that, first, the presence 
of A | qualitatively changes the dispersion relation of the 
transverse magnons, and second, the Eqs. I|3.5|l constitute 
an upper bound on this change, which becomes exact in 
the limit ^a/Cto ~* oo. We will discuss this point further 
in Sec. M below. 

Equations l|3.2l) for the longitudinal susceptibility in 
one-loop approximation remain valid. This expression, 
with xt,± from Eqs. (|3.5|l . and used in the Eliashberg 
equations, will provide an upper limit for the effect of 
A T on T cl . 



A. Linearized gap equation, and T c formula 



We start by linearizing the gap equations in Aj,j22 
while keeping the full dependence on Af. Equations 
become 



Ai(fc) 



r t / XL (g-fc)A i ( g )|G i (g)| i 



(4.1a) 



r t / 

Jq 

+2T t I [XT,+ (g-fc)+iXT,-(«-fc)] Gj 1 (q)/d 1 (q) , 

(4.1b) 

It / XL(q~k)A^(q)/d^q) , (4.1c) 



A T (fc) : 
E T (fc) : 

+2r t 



It / XL(q-k)Gj 1 (q)/d^q) 

Jq 

[XT,+ (q- k) +ixT,-(q- k)] G L (-q) , 



(4.1d) 



Notice that the down-spin and up-spin equations are cou- 
pled in two ways; directly via the explicit dependence of 
the transverse spin-fluctuation contribution to £j on Af , 
and indirectly via the dependence of the xt on Af . The 
former effect means that in the second contribution to 
Sf , the wave number is not strictly pinned to the up-spin 
Fermi surface. However, this effect is only on the order 
of Af /e F , and we neglect it. The up-spin and down-spin 
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equations then formally decouple, except for the depen- 
dence of x on Af. We can then solve for T^, in terms 
of Xi m the same McMillan-Allen-Dynes approximation 
that was employed in I for T c f. We find 



T ci = T (i)e-( 1+d °i +2d n)/ d w 



(4.2) 



Here To(t) is the same temperature scale that was used 
in I, viz., 6 



T (t) = T [e(t)t + S(-t)5\t\/4] 



(4.3) 



with To a microscopic temperature on the order of the 
Fermi temperature. The coupling constants d are given 
by 



d° 



dx x (1 — x 2 /2) XL^fepj,, iO) , 



dx a; xl (a^Fj. , iO) > 



v+ 



dx xxT,+(xk Fl ,iO) , 



(4.4a) 
(4.4b) 

(4.4c) 



with y± = kp^/kpi ±1 = S+/6- ± 1. Here ky a an d 
Nfo-, cr =t, |, are the Fermi wave number and the density 
of states at the up-spin and down-spin Fermi surface, 
respectively. 



B. The superconducting phase diagram 
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FIG. 6: T c j and T c | in zero-loop (dotted lines, right scale) and 
one-loop (solid lines, left scale) approximation. Parameter 
values used are 6l = 0.23, &t = 0.25. 



The calculation then proceeds as for T c p A representa- 
tive result is shown in Fig. We see that the relative 
difference between T c f and T c f in one-loop approxima- 
tion is on the same order as in zero-loop approximation, 
viz., at most about 10%. The theory thus predicts two 
superconducting transitions close to one another. In the 
notation of the schematic phase diagram in Fig. □ this 
means that the SCFM I phase is very narrow. 



For calculating T c numerically, we use the same mean- 
field relation between t and m or 5 as in I, namely^ 



t = 1 - (1 + 3r / 2 ) 1 / 3 /(l + ^? 2 /3) , (4-5a) 
m/n B n c = 3ry(l + 77 2 /3)/(l + 3r/ 2 ) , (4.5b) 

with n e the electron number density and /zb the Bohr 
magneton. In zero-loop approximation, with x given by 
Eqs. H3.2all3~5l) . T c j for a given A-f is now easy to calcu- 
late. Rather than solving the up-spin Eliashberg equa- 
tions for Af, we have approximated Aj ~ 2T c j. T c | is 
the same as in 1^ Our results are very similar to those 
obtained by Fay and Appel 6 and are shown in Fig. El The 
effect of Af on T c f is so small that it is not visible on the 
scale of the figure. At the zero-loop level, the feedback 
effect is thus unobservably small. 

At the one- loop level, with xl given by Eqs. (|3.2() 
with Eqs. (|3.5(l as input, the calculation is more com- 
plicated. We have approximated the effect of the mass 
in xt by a lower cutoff in the frequency summation, 
£}- = 2b T S{A 1 /S) 2 \n(S/T). This allows to perform the 
wave-number integral analytically, as was described in I. 
As in the zero- loop calculation, we replace Af by 2T c f. 
This procedure will overestimate the effect of Af on T c f . 



C. The specific heat 

The above results for the critical temperatures have 
implications for the specific heat. The two superconduct- 
ing transitions close to one another imply two features in 
the specific heat. At the mean-field level, the specific heat 
at a superconducting transition has a discontinuity, and 
for our choice of an order parameter the low-temperature 
specific heat is a quadratic function of the temperature^ 
One therefore qualitatively expects the specific heat coef- 
ficient as a function of temperature to behave as shown in 
Fig. Ufa). This should be compared to a situation where, 
as suggested in Ref.|3j, the down-spin electrons do not pair 
down to the lowest observable temperature, which would 
result in a specific heat coefficient as shown qualitatively 
in Fig. [7{b) . The data of Ref . show a very broad fea- 
ture, even for the best samples, which is consistent with 
either scenario. The prediction of the present theory is 
that, with further improving sample quality, the residual 
value of the specific heat coefficient will decrease, and the 
broad peak will be resolved into two discontinuities. 
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FIG. 7: The specific heat coefficient of URhGe as observed 
in Ref. (circies). The soiid line in (a) shows the qualita- 
tive theoretical prediction for a perfect sample based on the 
present analysis. The one in (b) would apply if the down-spin 
electrons remained unpaired to the lowest observable temper- 
ature. The size of the discontinuity is known theoretically 
only in the weak-coupling limitj— for the purpose of this qual- 
itative illustration we have adjusted it to fit the maximum in 
the data. 



V. DISCUSSION AND CONCLUSION 

Let us further discuss several aspects of our theory and 
its implications. First, we consider in more detail our re- 
sults regarding the magnetic susceptibility. In a normal 
conducting ferromagnetic phase, the transverse spin sus- 
ceptibility is soft, displaying ferromagnetic magnons with 
a quadratic dispersion relation. The scale in the disper- 
sion relation is set by the Stoner gap or exchange splitting 
5, which is proportional to the magnetization. In natural 
units, measuring wave numbers in units of the magnetic 
coherence length £ m , which is on the order of the Fermi 
wave number fcp, we have 

At zero frequency, the dimensionless transverse magnetic 
susceptibility diverges for fc — > like 

XT (fe) =c/(kU) 2 ■ 

For a given distance from the magnetic phase transition, 
c is a number of order one whose exact value depends 
on the details of the model. The Gaussian approxima- 
tion we have employed correctly reflects these exact fea- 
tures, see Eqs. Q3.1|l . In the presence of an up-spin gap, 
the structure of the magnon changes qualitatively. The 
Gaussian approximation for xt predicts a true gap, Eqs. 
(|3.5|) . This is a consequence of the fact that the Gaussian 
approximation effectively treats the superconducting gap 
as a fixed external field. The LGW theory of Appendix 
shows that this does not correctly reflect the symmetry 
properties of the system, and that the static \t actually 



has the structure 



XT(fe) 



1 



9- 



2A F r t |A T | 2 /5 2 



(5.1a) 



Here g is related to the coupling constant g and the mag- 
netic order parameter rh of the LGW theory by g = g rh. 
This result shows that, in general, xt has a much more 
complicated structure than in the absence of spin-triplet 
superconductivity. At asymptotically small wave num- 
bers, fc 2 < <?£a 2 ' g drops out of the expression for xt, 



and we have 



fc 2 l + {^/5f{UIU? 



(5.1b) 



Restoring the frequency, the dispersion relation of the 
magnon then is 

[l + (A T /<5) 2 (£ A /£ m ) 2 ] fc 2 . 

We see that the magnon dispersion relation is substan- 
tially changed by the presence of Af as long as the con- 
dition expressed by Eq. (|A10|) is fulfilled. Namely, the 
magnon is much stiffer, i.e., the frequency rises much 
faster with the wave number, than in a normal conduct- 
ing ferromagnet. The magnitude of this effect will depend 
on the detailed parameters of the material in question, 
but it clearly can be substantial: Assuming a value of Aj 
on the order of the superconducting T c , or Af w IK, a 
value of 8 on the order of 10 times the magnetic T c j24 or 
S m 100K, and £a/£™ ~ 1,000, we find that the pref- 
actor in the dispersion relation is enhance by a factor of 
100 over its value in the normal conducting phase. 



For larger wave numbers, k > g£ A , we have 
Xt(*0 



(kU) 2 + g2N F r t \Atf/6 2 



(5.1c) 



Using Eqs. I|A9J| . we have the correspondence g = 
g(S/2T)/3 > 1, so we recover the result of the Gaus- 
sian theory, Eqs. (|3.5() . As long as Eq. (jAlOfl is valid, 
Xt thus displays a pronounced plateau as a function of 
the wave number, followed by a very steep increase at 
asymptotically small wave numbers. The Gaussian the- 
ory approximates this behavior by a true gap, ignoring 
the asymptotic regime. In the limit £A/£m oo the 
asymptotic regime shrinks to zero, and the Gaussian the- 
ory becomes qualitatively correct. However, in this con- 
text one should note that the LGW theory of Appendix 
lAlis valid only for wave numbers |fe| < £7 . Within the 
LGW theory, the shoulder will lie in that regime provided 
the coupling constant g is sufficiently small. While the 
above identification of g with the parameters of the mi- 
croscopic theory makes g effectively of order unity, given 
the initial sharp rise of the magnon frequency with the 
wave number in the regime where the LGW is valid, a 
pronounced shoulder in the dispersion relation at inter- 
mediate wave numbers is inevitable. The wave number 
region where one expects this shoulder is given by 



9&<k 2 



<<7(|A T |/<5) 2 & 



(5.2) 
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If we use the same numbers as above, and g = 1, £ m = 
lA, the upper limit of this wave-number range is given by 
|fe| = O.OlA. This is a factor of 3 below the smallest wave 
numbers currently observable with neutron scattering. 25 
However, in materials with smaller values of the exchange 
splitting 8 the plateau should be in an observable regime. 

Second, we add some comments about our prediction of 
two superconducting transitions that are separated only 
by a small temperature interval. Within our model, we 
have found this prediction to be very robust, especially 
given that all of our approximations have a tendency 
to overestimate the suppression of T c j. compared to T c |. 
Even an (artificial) increase of the gap in the transverse 
susceptibility by a factor of 10 does not make a visible 
change in Fig. This reflects the fact that the up-spin 
and down-spin pairing are mediated by the same effective 
potential, viz., \l, see Eq. (|4.4a|l and the corresponding 
Eq. (3.14a) in I. While \l is modified by Af, the effect 
is not sufficiently large to lead to a substantial separa- 
tion of T c | and T c j. This in turn means that the lower 
value of T c i is overwhelmingly due to the lower value of 
the density of states at the down-spin Fermi level. In 
this context we need to keep in mind that we have used a 
free-electron model with parabolic bands. A complicated 
band structure could lead to a drastically reduced value 
of iVpj, which in turn would lead to a much lower value 
of T c i. If experiments on samples of improved quality 
should fail to show two transitions, this would be the 
most likely explanation. This would be of interest also 
with regard to distinguishing between the two proposed 
explanations for why the superconductivity is observed 
in the ferromagnetic phase only: Sandeman et ali2& have 
proposed a mechanism based on an intricate structure of 
the density of states, while the explanation proposed in 
I is based on properties of the magnetic susceptibility. 

Third, we come back to our discussion of the specific 
heat. The experiment of Ref. shows that the observed 
superconductivity is indeed a bulk effect. What is not 
clear a priori is the origin of the large residual value of 
the specific heat coefficient. While the down-spin elec- 
trons remaining unpaired, as was suggested in Ref. 0, is 
a possibility, normal-conducting regions within the sam- 
ple would have the same effect and would also lead to 
the observed smearing of the discontinuity in the spe- 
cific heat. With increasing sample quality, the discon- 
tinuity should become sharper, and a crucial question 
will be whether the residual value drops correspondingly. 
Of course, the emergence of the predicted double fea- 
ture from the narrowing peak will be the most direct 
test of our predictions regarding the specific heat. In 
this context it is interesting to note that such a split 
transition, with two closely spaced discontinuities in the 
specific heat, has been observed in UPtajSL 2 ^ but only af- 
ter a long period of gradually increasing sample quality. 
(UPt3 is not ferromagnetic, though, and the split tran- 
sition has a physical origin that is very different from 
what we have discussed.) Even in the best samples that 
show the split transition, however, there is a substantial 



residual specific heat coefficient, the origin of which is 
not quite clear. In a ferromagnetic superconductor, one 
also has to keep in mind that the ground state will not 
be homogeneous, due to the formation of a spontaneous 
vortex stateiS2i22i2i Normal electrons in the vortex cores 
are a possible source of a residual specific heat coefficient. 
This effect has been neglected in the current theory and 
will be pursued in a future publication. 

Fourth, we briefly comment on the fact that the mass, 
or pseudo-mass, induced in the transverse magnetic sus- 
ceptibility by the superconductivity, is a singular func- 
tion of the magnetization. This can be seen in Eqs. (|3.5|l . 
and it also leads to the factor of 1/6 in the relation be- 
tween the phcnomcnological coupling constant g in the 
LGW theory and the microscopic parameters, Eq. (|A9c|l . 
This behavior can be traced back to the behavior of the 
integral in Eq. I|3.5c[) . and thus ultimately to the Green 
functions and the soft particle-hole excitations that are 
characteristic for itinerant electron systems. The sin- 
gularity is therefore a result of a coupling between the 
particle-hole excitations and the magnetic and supercon- 
ducting Goldstone modes. It is very similar in nature 
to, e.g., the anomalous magnetization dependence of the 
magnon stiffness in a normal conducting ferromagnetic 
phase that was discussed in Ref. 133 . 

In conclusion, we summarize our results. We have 
presented a consistent and self-contained theory for the 
coexistence of superconductivity and ferromagnetism in 
itinerant electron systems. We have presented a field- 
theoretic formulation of this problem that allows for the 
determination of both the magnetic and the supercon- 
ducting equation of state in a systematic loop expansion. 
This method, which utilizes Ma's procedure for gener- 
ating equations of state, remedies some shortcomings of 
the earlier theory presented in I, which relied on a mini- 
mization of the free energy The self-contained character 
of the theory is achieved by means of explicit expressions 
for the magnetic susceptibility, which is needed as input 
for the generalized Eliashberg equations. These expres- 
sions have been evaluated in ferromagnetic phases, both 
normal conducting and superconducting ones. This is a 
generalization of the theory for the magnetic susceptibil- 
ity that was developed in I, and it explicitly takes into 
account the feedback effects that are characteristic for 
any purely electronic mechanism for superconductivity. 

We have explicitly evaluated this theory to one-loop 
order, and for a model that allows for two components 
of the superconducting order parameter, one each for 
Cooper pairs consisting of up-spin electrons and down- 
spin electrons, respectively. The limitations of this 
one-loop approximation, which neglects superconduct- 
ing fluctuations and uses a zero-loop expression for the 
magnetic susceptibility, have been discussed by means of 
a phcnomcnological LGW theory that complements our 
microscopic theory. We have found that, for generic pa- 
rameter values, the two superconducting transitions that 
describe the pairing of up-spin and down-spin electrons, 
respectively, occur close to one another, with transition 
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temperatures that typically differ by only on the order 
of 10%. This suggests that, if the superconductivity ob- 
served in URhGe is indeed of a spin-triplet p-wave type 
mediated by ferromagnetic fluctuations, then the broad 
feature observed in the specific heat near the temperature 
of the resistive transition^ should contain two transitions 
that are close together. If samples of improved quality 
should show only one sharp discontinuity in the specific 
heat, then this would be a strong argument against the 
type of pairing we have assumed in this paper. A caveat 
is provided by our assumption of parabolic bands, how- 
ever, as discussed above. The presence of spin-triplet 
superconductivity has further been shown to drastically 
change the structure of the dispersion relation of the fer- 
romagnetic magnons. In materials with a small exchange 
splitting this effect is observable with neutron scattering, 
and can also be used as a probe for the nature of the 
superconductivity. 
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APPENDIX A: LGW THEORY 

In this appendix we discuss the mean- field action A^ ' , 
Eq. 1)2. 8[l . or the corresponding Landau free energy den- 
sity f(°> = — (T /V)A^°\ and its generalization to a phe- 
nomenological LGW theory. 

1. Landau theory 

From Eq. we have 

Z (0) = y™ 2 -f TrtA^-^TrlnG- 1 . (Al) 

In what follows we neglect the normal self-energy con- 
tribution to the matrix A, Eqs. I|2.6|l . which lead only to 
a trivial renormalization of the normal Green function. 
An expansion in powers of the order parameters m and 
A then yields 

/(0) = Tr^O-Go/^-^Tr^Go) 2 

+ y 1 ™ 2 ~ ~w m Tr r° 73 <5oAT 60 aT ) + ■ ■ ■ (A2) 



Here t = 1 — 2N-pT t is the mean-field distance from the 
magnetic critical point, and we have omitted cubic terms 
of order to 3 , to 2 A, and A 3 , as well as all quartic terms. 

This mean-field free energy does not describe a super- 
conducting transition; this requires magnetic loops, as 
we have seen in the main text. The magnetic part, how- 
ever, agrees with an expansion of the mean- field magnetic 
equation of state, Eq. (|2.14|) . In particular, the cubic 
term of 0(toA 2 ) shown in Eq. (|A2|) has important conse- 
quences. For a nonzero superconducting order parameter 
A T , Eq. p.fijl . the quantity 



7 1 / \ 

h eS = ^ Tr (£o 73 G A T G A T J (A3a) 

acts like an effective external magnetic field. Performing 
the trace, we have 



/ l e ff =-r t -^|A T (fc)| 2 |Go(fc)| 2 G (fc) . (A3b) 

fc 

General arguments^ show that such a term in the free 
energy leads to a mass in the magnetic Goldstone mode, 
i.e., in the transverse spin susceptibility, 

X - 1 (k = 0) = h cS /m . (A4) 

If we expand Eq. (|3.4|) in powers of Aj and to, we see 
that this is indeed the same result we obtained in Sec. 
IIIII by a direct calculation of the susceptibility. The cur- 
rent derivation makes it obvious that this result has been 
obtained while neglecting fluctuations. We consider the 
influence of fluctuations in the next subsection. 



2. LGW theory 

The conclusion in the previous subsection, namely, 
that the transverse spin susceptibility is massive, cannot 
be strictly correct. Even in the presence of superconduct- 
ing order, one expects the system to be invariant under 
rotations of all spins, and this symmetry must manifests 
itself via a Goldstone mode to which the magnetic sus- 
ceptibility must couple. To investigate this point, we 
consider the following phenomenological action^ 
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S[m, 4>] = I dx t m m 2 (x) + ^ n (Vm(x)) +u m m A (x) 
-ig I dx m(x) ■ (4>(x) x <p*{x)) 



dx [t A \<f>(x)\ 2 + e A \V<I>(X)\ 2 + «A \<t>{x)\ 4 



(A5) 



Here m and <p are the magnetic and the superconduct- 
ing order parameter, respectively, and the latter has been 
represented as a complex 3-vector in spin spaced t m and 
t a are the dimensionless distances from the magnetic and 
superconducting critical point, and £ m and £a are the 
(zero-temperature) magnetic and superconducting coher- 
ence lengths, respectively, g is a phenomenological cou- 
pling constant. At the mean-field level, m = (0,0, m), 
and 4> = (A,iA,0), with rh oc m and A c>c Aj, and the 
theory is seen to have the structure of the Landau theory 
derived from the miscroscopic theory in Sec. IA ll above. 

This action is invariant under a symmetry group 
SO(3)xU(l), which represents co-rotations of the vec- 
tors m and <p, and in addition a gauge transformation 
of the complex field 0. In the ordered state given above, 
this symmetry is spontaneously broken, and the ordered 
state is invariant only under rotations about the 3-axis. 
We thus expect dim(SO(3)xU(l)/SO(2)) = 3 Goldstone 
modes. An explicit calculation, involving an expansion to 
Gaussian order about the ordered state, confirms this ex- 
pectation. Of the three soft modes, two are spin wave-like 
modes, namely, superpositions of the transverse compo- 
nents of m and the 3-component of <p, while the third 
is an Anderson-Bogoliubov-like mode that reflects the 
broken gauge symmetry. Specifically, the transverse m- 
susceptibility has an overlap with one of the spin-wave 
Goldstone modes. The Gaussian theory yields 



and the function g(S/2T) in Eq. i|3~5c|) . 



Xt(p) = (mi(p)mi(-p)) 



(5^ A 2 +P 2 ) 



2p 2 s(A^ 2 + m^-2 



rhp 2 



(A6) 



From this result we see that, in the momentum range 

jmG 2 <P 2 <9(A 2 /«)C 2 , (A7) 
the transverse susceptibility is effectively massive, 

X T (p)«m/2ffA 2 . (A8) 

m and A can be related to m and Aj, respectively, 
in the microscopic theory by comparing Eq. I|A5(1 with 
the Landau theory in Sec. IA II above. One finds the 
correspondence 



m = m^T t /2T , 
A = A T v/Af/T . 



(A9a) 
(A9b) 



g does not corresponds to a simple constant in the micro- 
scopic theory, but rather to a function of m (or S = T t m) 
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^2TT t — g(S/2T) 
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(A9c) 



This complicated behavior of the effective coupling be- 
tween the magnetic and superconducting order parame- 
ters is not reflected by the LGW theory. It is a result 
of the itinerant nature of the electrons, and the corre- 
sponding soft particle-hole excitations, as is discussed in 
Sec. |V| Upon substituting Eqs. (|A^|) in Eq. (|A"8|) we see 
that the latter is the same result as Eq. I|3.5a|l for the 
magnetic susceptibility. 

The momentum range given by Eq. (|A7|I exists only 
if the condition A ^ m£, m /£,A is fulfilled. In terms of 
the Stoner gap S = T t m and Af this condition takes the 
form 

£/A T <W£m ■ (A10) 
These results are further discussed in Sec. 



APPENDIX B: THE SPIN SUSCEPTIBILITY 

In this appendix we establish the relation between the 
physical spin susceptibility \s, as measured, e.g., by neu- 
tron scattering, and the quantity \ that emerges as the 
SM propagator of our field theory, see Eqs. p. 91 12.l0|l . 
We start by adding a magnetic field h(x) to our action 
that couples linearly to the electron spin density and acts 
as a source field for spin-density correlation functions. 
The action, Eq. (|2.1I) . becomes 

A[M,g, k;h] = Tr (AG) - J dx M(x) ■ M (x) 
+ iTrln(G 1 + ■ M - A T + ~f ■ h) . (Bl) 

The partition function, 

Z[h}= J D[SM,5\,6G] e A[M ' s ' k ' h] , (B2) 

then serves as the generating functional for the magneti- 
zation, 



m = (n^(x)) 



Sh 3 (x) 



In Z[h] 



(B3) 



h=0 
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and for the spin susceptibility, 

Xs(x-y) = (Sn % a (x) 8n J s (y)) = 



5hi(x) Shj(y) 



\nZ[h] , 



h=0 



(B4) 



where n s (x) is the electron spin density, and Snl(x) = 
K(x) - (ni(x)). 

The simplest way to deal with the source held is to 
shift M: M(x) -> M(x) - h{x)/y/W~ t . We then find 



m=y/2/Tt{M z (x)) 



(B5) 

in agreement with Eq. (|2.5a|) . An evaluation of (M 3 (x)) 
to one-loop order yields the magnetic equation of state, 
Eq. I|2.14[) . Similarly, 

X« (x-y)= (SMi(x) 8M j (y)) ~ % S(x - y) . (B6) 

If we evaluate the 8M -correlation function in Gaussian 



approximation, we fim 
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Xl j (x - y) w [xy(a; - y) - % <5(x - y)] /r t 



(B7) 



We see that the ^TW-propagator of the field theory is 
simply related to the physical spin susceptibility. In par- 
ticular, the transverse part of Xs is massive if and only if 
that of x is. For the Fourier transform of the longitudinal 
part, we find from Eq. 



X^(k) = X ?(k) 



-xl(k) 



l + TtxUk) 



(B8) 



which has the expected RPA-type structure. Notice that 
the "contact term", which results from the term in A 
that is quadratic in h, provides the numerator in Eq. 
(|B8|) that is missing in xl- 
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